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ABSTRACT 

We are undertaking a large scale radial velocity survey of the Galactic bulge which uses M giant 
stars selected from the 2MASS catalog as targets for the CTIO 4m Hydra multi-object spectrograph. 
The aim of this survey is to test dynamical models of the bulge and to quantify the importance, if 
any, of cold stellar streams in the bulge and its vicinity. Here we report on the kinematics of a strip 
of fields at —10° < / < +10° and b = —4°. We construct a longitude-velocity plot for the bulge 
stars and the model data, and find that contrary to previous studies, the bulge does not rotate as a 
solid body. From —5° < I < +5° the rotation curve has a slope of roughly 100 km s _1 kpc -1 and 
flattens considerably at greater / and reaches a maximum rotation of 45 km s -1 . We comp are our 
rotation curve and velocity dispersion profile to both the self-consistent model of (|Zhaolll99a ) and to 
N-body models; neither fits both our observed rotation curve and velocity dispersion profile. The high 
precision of our radial velocities (~ 3 km s _1 ) yields an unexpected result: hints of cold kinematic 
features are seen in a number of the line of sight velocity distributions. 

Subject headings: Galaxy: bulge - Galaxy: kinematics and dynamics - Stars: late-type - 
stars:kinematics - techniques: radial velocities 



1. INTRODUCTION 



The COBE 2/zm image of the b ulge (|Dwek et al.lll995t 
iBinnev. Gerhard, fc Spe rgel 1997]) and models of the pro- 
jected 2/Ltm light shows a bar-like structur e that is also 
detec ted in star counts of red clump stars (jStanek et al.l 
119971 ). The anom alously high optical depth of the bulge 
to microlensing (|Alcock et al.l [2000) can be explained 
only by assuming a bar whose major axis extends roughly 
toward the Sun, thu s raising the rate of star-star events 
(|Han fc Gouldl 120031 ). Theoretical mod els of the b ulge 
initially followed axisymmetric models (jKentJ IT992T ) but 
have g raduated to self-consistent rapidly rotating bar s 
(|Zhad Il996t lHafner et all 120001 : iBissantz et all l2004h . 
with the density and the potential strongly constrained 
by the observed microlensing rates in the bulge as well as 
gas kinematics. However, the phase space of the bar is 
relatively incompletely constrained by stellar kinematic 
data. 

Study of the kinematics of the bulge is complicated by 
the large and variable foreground extinction, the pres- 
ence of a contaminating disk population extending from 
the foreground well into the Galactic Center, and source 
confusion arising from the high density of stars. Posi- 
tional measurements from wide field Schmidt plates are 
therefore impossible, where source confusion makes any 
astrometric exercise daunting. 

M giants, while faint in the traditional optical band- 
passes due to their cool temperatures and TiO bands, 
are bright in the / band (13 < I < 11) and are easy tar- 
gets for spectroscopy, if positions are known. They fur- 
ther have the advantage of being ubiquitous throughout 
the bulge and are luminous enough to be studied even 

1 Department of Physics and Astronomy, UCLA, Los Angeles, 
CA 90095-1547 rmr@astro.ucla.edu 

2 SUPA, School of Physics and Astronomy, University of St An- 
drews, KY16 9SS, UK 

3 Visiting Astronomer, Cerro Tololo Inter- American Observa- 
tory. CTIO is operated by AURA, Inc. under contract to the 
National Science Foundation. 



in fields with substantial extinction. Finally, the short 
lifetimes of AGB stars and luminous giants limit their 
n umbers e nough that source confusion is not an issue. 

iMouldl dlflSah was the first to measure the veloc- 
ity dispersion of the bulg e in Baade's Window, based 
on the M giant sample of Bl anco. McCarthy, fc Blancol 
(p8l that was first classified from low dispersion slit- 
less spectra. Multi-fiber spectr o scopy of this sam- 
ple ([Sharpies. Walker, fc Cropper! [T990) subsequently 
harvested roughly 250 bulge giant velocities. De- 
spite numerous investigations of the dyna mics of 
stars in the dire c tion of Baade's Window (fMouldl 
1981 [Richl H990l ISharples. Walker, fc Cropped 11995 
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dTvson fc Rich 19 9 H iMinniti et alJ Il992t iMinnitil 119961 : 
iBlum et al.l [l99 4. 1995) there has been, up to now, no 
large scale survey of the dynamics of th e stellar popula- 
tion in the bulge. iBeaulieu et al.l ([20001 ) survey the plan- 
etary nebula population; both the rarity of PNe (due 
to their brief lifetimes), the problem of disk or bulge 
membership, and the considerable distance uncertainty 
make PNe a problematic population; We believe a well 
selected sample of M giants is likely the best probe of the 
b ulge/bar population . 

iFrogel fc Whitfordl (|1987f ) showed that the M giant 
luminosity function is consistent with an old popula- 
tion roughly the age of the halo and globular clus- 
ters. The first detaile d abundance analysis of M giants 
(jRich fc Origlial [20051 ) finds a n abundance range simi- 
lar to that of the K giants (|Mc William fc Richl [l99l 
iFulbright et al.|[2~006l ). In principle, the most metal-poor 
component of the bulge population might not evolve 
through the M giant phase as the K giant abundance 
distribution is peaked at slightly sub-solar metallicities, 
but most of the stars are more metal-rich than 47 Tuc 
(at —0.7 dex) and are therefore candidates to reach the 
luminosities a nd effective temper atures characteristic of 
the M giants (jZoccali et al.ll2003fl. 

The proper motion study of ISumi et all (|2004l ) ad- 
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Fig. 1. — Fig. 1: Color-magnitude diagram of 2MASS candidates 
and filled symbols (o bserved stars), including reddened isochrones 
HGirardi et al.l 12003) for [F/eH]=-1.3 and -0.5. The parallelo- 
gram indicates our selection region; the blue cutoff rejects many 
objects that are closer than the bulge, which have lower reddening 
and are brighter than the red giant branch. The reddening vector 
corresponds to E(J - K) = 0.33 from the ISchlegel et~aLl llT998h 
map. 

dresses a number of fields in the bulge (avoiding high 
extinction) and eventually, a second epoch of 2MASS 
imaging would provide proper motions for large numbers 
of giants. However, the addition of radial velocities and 
ultimately metal abundances, is needed for a complete 
dynamical model and the M giant population provides 
the perfect sample of stars to target. 

The dynamical model for the bulge/bar has a number 
of important implications. Large samples of uniform ra- 
dial velocity data are still of great value in constraining 
the bar vs. axisymmetric models, and the nature of the 
orbit families supporting the bar. Further, the interpre- 
tation of the microlensing events in the bul ge depends on 
the u se of an accurate dynamical model (|Han fc Gouldl 
2003). The recent discovery of pla netary transit host 
stars in the bulge (|Sahu et al.i r2006) gives an additional 
incentive to improve our knowledge of the bulge/bar 
model. 

With the availability of the Two Micron All-Sky Sur- 
vey 2MASS (jSkrutskie et al.ll2006l ) we realized that the 
key ingredients of high precision positions and photome- 
try would finally be available everywhere in the Galactic 
bulge except near the plane of the Galaxy. At this time, 
we present ~ 2300 spectra and have obtained a radial ve- 
locity precision of ps 3 km s _1 for our most recent (2006) 
data. 

In the past, optical radial velocity studies in the bulge 
have not emphasized high precision, because of the large 
velocity dispersion and the expectation that the short 
orbital periods would phase-wrap any cold structures 
out of existence in well under a Gyr. Yet our preci- 
sion is sufficient to enable a search for more cold streams 
analagous to those associated with the Sagittarious dwarf 
spheroidal galaxy; some candidate cold features are seen 
and followup observations are underway (Reitzel et al. 
2007 in prep.). Here we report on the kinematics of stars 
along a band at b = —4°, obtaining a rotation curve and 
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Fig. 2. — Fig. 2a-c: Velocity dispersion profile of three minor 
axis fields ranging from —3° < b < —4°; note the dumpiness in a 
and b. The velocity dispersion profile of the sum of these fi elds is 
given in Fig 2d. The dotted line is the prediction from the IZhaol 
( 1996) model. 

velocity dispersion profile. 

2. OBSERVATIONS AND SAMPLE SELECTION 

The choice for optical spectroscopy is driven by the 
availability of multi-object wide field sp e ctrosc opy. IR 
spectroscopy, such as that of iBlum et al.l (|1995h . can be 
used to great effect in the fields of highest reddening but 
obviously yields far fewer spectra. Here we present a brief 
sketch of the selection method and analysis; full details 
will be given in Howard et al.(2007 in prep). 

We use the sur vey method of 

ISharples. Walker. fc Cropperl (|1990[ ) as the model 
for our program. They observe that stars with I < 11.8 
have a lower velocity dispersion and are likely disk 
members. When the same field is examined in 2MASS, 
the K vs J — K CMD shows a clearly defined red 
giant branch. In the SCW90 field, we find that the 
magnitude limit of / < 11.8 corresponds roughly to 
K < 8. Further, there is break in the luminosit y 
function at this K magnitude (jFrogel fe W hitford 1987). 
While a few bulge AGB members may be present at 
brighter magnitudes, the bulk of such bright stars will 
be foreground contaminants. We adopt a range of 
9.25 < K < 8 as the basis for selecting the 2MASS M 
giants (Figure 1). 

A parallelogram-shaped selection region is adjusted in 
color and magnitude to encompass completely the red 
giant branch locus of old stars at the distance and red- 
de ning of the fie l d as indicated by the reddening map 
of ISchlegel et all (|1998t ). We find that reddening varies 
greatly, ranging from 1.5 < Ay < 5.0. In fields with high 
extinction, differential reddening is an issue and we widen 
the selection region to accommodate the red giant branch 
as observed (Figure 1). The indicated isochrones show 
that we admit stars with —1.3 <[Fe/H]< 0.3 a range that 
spans the entirety of the lZoccali et al.l (|2003l ) abundance 
distribution 

2.1. Spectroscopy 

We use the Hydra multi-fiber spectrograph at the 
cassegrain focus of the Blanco 4m telescope at Cerro 
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Fig. 3. — Fig. 3a: Rotation curve fro m our data (filled symb ols); 
open symbols indicate binned PNc from Bcaulieu ct al. (2000); the 
PNe agree reasonably well with our data and the departure from 
solid body r otation is clear. The solid line in dicates the model of 
IZhad {1996) while the dotted line is that of IFuxl (119971 ) and the 
dashed line is from Scllwood (1993) (sec text). Fig. 3b: Velocity 
dispersion is indicated with the PNe and models as above. 

Tololo. We optimize for spectroscopy in the red, tak- 
ing advantage of the red colors of M giants and employ 
the KPGLD grating, blazed at 8500A giving 0.45A pix -1 
with 2-pixel on-chip binning yielding an effective resolu- 
tion of 0.88A pix -1 and a full spectra range of 1800A. 
We use the 200/Ltm slit plate, giving us an increase in 
resolution; the loss of light from the slit plate is incon- 
sequential to our S/N because the stars are so bright. 
Our useful spectral range in 2006 was 6891A to 8714A; 
the first two lines of the Call infrared triplet are often 
well detected and contribute to our success in obtaining 
a correlation peak. However, there are also stars that 
are so red that the Ca triplet is overwhelmed by TiO 
absorption; in these cases, the wide wavelength range is 
essential for the velocity measurement. Our exposures 
are typically 3 x 900 s but longer exposures were needed 
when cirrus was present. Each field has on average, 108 
successfully exposed M giants and 20 fibers that are used 
to obtain a sky background spectrum. Flat fielding, sky 
subtraction, throughput correction, scattered light sub- 
traction, and wavelength calibration for each exposure 
were all accomplished using the IRAF task dohydra. The 
spectra are binned to 34.3 km s _1 pix -1 and normalized. 
Radial velocities are measured using the fxcor task in 
IRAF; this requires the spectra to be Fourier filtered ex- 
clude features greater than 50 pixels or smaller than 3 
pixels in extent. Regions of the spectrum with obvious 
telluric features were excluded, leaving only 60% of the 
spectrum usable for cross correlation. Our final cross 
correlation regions are 7000-7150A, 7300-7580A, 7700- 
8100A, and 8300-8600A. 

3 standard stars are all used in the cross correlation, 
with the final velocity for any given star being the av- 
erage of the 3 derived velocities. The three standards 
return velocities that agree to within 1.6 km s _1 on av- 
erage, with a standard deviation of 1.4 km s _1 in these 
differences. 
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Fig. 4. — Fig. 4: Upper panel: longitude-velocity diagram for 
the data along b = —4°, which have been smoothed to 1° bins and 
by 8 km s — 1 . Lower Panel: The same region extracted from the 
IZhaol 1119961 ) model. The lower right region would be populated by 
stochastic orbits in the model. By construction, the model has no 
retrograde orbits; their inclusion would help bring the model into 
agreement with the data. 

We report radial velocities for a total of 2294 M giants. 
Here we consider the fields spanning across b = —4°. 
We compare our observed velocities and velocity disper- 
sion with those pre dicted by th e self consistent rotating 
bulge/bar model of IZhaol (| 19961 ). 

3. THE ROTATION CURVE AND VELOCITY DISPERSION 
PROFILE 

We now discuss the results from our study, beginning 
with the minor axis velocity dispersion (Figure 2) re- 
sulting from the sum of all fields on the minor axis at 
b = —2.5° to —6.5° as well as three examples of the 
fields that went into the summed distribution. The best 
fit Gaussian to the data gives a = 119 ± 5 km s _1 with 
vq = — 11 ± 30 km s . A number of peaks are present in 
this histogram; these are prominent in the contributing 
histograms (Figure 2a-c) . A c ompar i son w ith the equiva- 
lent region extracted from the I Zhao! (|1996l ) model is also 
shown; the agreement is very good, but the velocities in 
the data appear to be more clumpy than they are in the 
model. More striking peaks (w 2.5a) are are found in 
some of our other fields, but because such features can 
occur in random draws, new observations are required 
for confirmation (Reitzel et al. 2007). 

Figure 3a, b shows our rotation curve and velocit y dis- 
persio n profile compared with that predicted by IZhaol 
(|1996h : we give the data in Table 1. We do not confirm 
the solid body rotation claimed by a number of previous 
studies; after reaching an amplitude of ~ 40 km s _1 , the 
rotation curve flattens beyond \l\ > 3°. The velocity dis- 
persion profile remains higher than 75 km s" 1 even for 
the fields at I = 10°. We c ompare our data with the plan- 
etary nebulae (PNe)from (jBeaulieu et alfe OOOh In order 
to have a reasonable number of PNe, we have accepted 
those in the range —8° < b < —3° and we have binned 
the data. Considering the less secure distances and as- 
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signment of population for the PNe, the agreement is 
good, and settle s the questio n of solid body rotation for 
the bulge. The IZhaol (| 19961 ) model is a self-consistent 
rapidly rotating bar that is constrained to have no ret- 
rograd e orb its. Also plotted are the N-body bars of iFuxl 
(|1997t ) and lSellwoodT(|1993l ). both of which are N-body 
bars formed from initially unstable disks. Sellwood's bar 
starts from a rigid Plummer sphere with a live Kuzmin 
disk, with mass ratio 3:7 and no dark matter. Fux's bar 
starts from an equilibrium of dark halo, power law nu- 
cleus, and an exponential disk; in contrast to the Zhao 
model that is fit to the data, a best-fit model is selected 
from a large number of N-body realizations. 

We adopt t he bar angles and ma ss normalisations as 
suggested by iBeaulieu et al.l (|2000f ) that give an over- 
all good match with the appearance of the COBE map. 
We project these models and calculate the velocity mo- 
ments in each line of sight without distinguishing disk 
and bulge particles. Neither model is satisfactory: Fux's 
best model gives a remarkably good fit to the rotation 
curve but has too low V ro t/cr. Sellwood's bar appears 
to have a very high V ro t/cr. It appears that our data is 
challenging to fit by both disk-instability formed bars as 
well as Schwarzschild models wit h fixed poten tials. 

In Figure 4, we compare the IZhaol (|1996l ) model in 
longitude-velocity space with our radial velocities, both 
sampling the slice at b = —4°; we have smoothed the 
data in longitude by 1.0° and 8 km s _1 in velocity. The 
addition of retrograde orbits may improve the agreement 
of the model and data; we are working toward this end 
(Zhao et al. 2007 in prep.). 

While our claim of slower rotation for the bulge ap- 
pears to contradict earlier studies, we emphasize our 



agreement with the PNe. It is also interesting that 
th e kinematics of K g iants in the (/,&) = (8,7) field 
of iMinniti et alj (|1992t ) agrees with our rotation curve 
(44.5 km s _1 ); faster rotation (and a lower velocity dis- 
persion) is observed in his (I, b) = (12, 3) but th at field 
that might be disk dominated. iMenziesI (| 19901 ) used a 
small number of Miras to sample velocities over a wide 
range in longitude, and found rapid solid body rotation. 
Small numbers, or contamination by members of the 
disk, might be invoked to explain the more rapi d rotation 
observed for SiO masers (jlzumiura et al.lll995[ ). 

Our data are clearly at odds with the widely held view 
that the bulge rotates as a solid body. The relatively 
high velocity dispersions of our fields is reassuring in the 
sense that the M giant selection criterion is yielding good 
kinematic probes of the bulge/bar population. Extend- 
ing this survey to a larger number of fields in the bulge 
offers the possibility of undertaking detailed tests of dy- 
namical models, something that is not presently possible 
in distant galaxy populations. This will give new insights 
into the structure and dynamics of the bulge, and into 
the formation of the Milky Way. 
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TABLE 1 

Observed Bulge Rotation Curve and Velocity Dispersion Profile 
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